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ABSTRACT 

In spatially distributed multiuser antenna systems, the received 
signal contains multiple carrier-frequency offsets (CFOs) arising 
from mismatch between the oscillators of transmitters and receivers. 
This results in a time-varying rotation of the data constellation, 
which needs to be compensated at the receiver before symbol re- 
covery. In this paper, a new approach for blind CFO estimation 
and symbol recovery is proposed. The received base-band signal 
is over-sampled, and its polyphase components are used to formu- 
late a virtual Multiple-Input Multiple-Output (MIMO) problem. By 
applying blind MIMO system estimation techniques, the system re- 
sponse can be estimated and decoupled versions of the user symbols 
can be recovered, each one of which contains a distinct CFO. By 
applying a decision feedback Phase Lock Loop (PLL), the CFO can 
be mitigated and the transmitted symbols can be recovered. The esti- 
mated MIMO system response provides information about the CFOs 
that can be used to initialize the PLL, speed up its convergence, and 
avoid ambiguities usually linked with PLL. 

^:e>'woriii'-Multi-user Systems, Distributed Antenna Systems, 
Carrier Frequency Offset, Blind MIMO System Identification 

I. BACKGROUND 

In both wireless and wireline communication systems, received sig- 
nals are often corrupted by carrier-frequency offsets (CFOs), due 
to Doppler shift and/or local oscillator drift. The CFO causes a 
time-varying rotation of the data symbols, and thus before symbol 
recovery, it must be estimated and accurately compensated for by 
the receiver. The CFO can be estimated via the use of pilots sym- 
bols; however, even a small error in this estimation tends to cause 
large data recovery errors. This necessitates transmission of pilot 
symbols rather often. In single user systems, or in multiuser sys- 
tems where the transmitters are physically connected to the same 
oscillator, there is only one CFO that needs to be estimated. This 
is typically done via a decision feedback Phase Lock Loop (PLL) 
at the receiver. The PLL is a closed-loop feedback control system 
that can adaptively track both frequency and phase offsets between 
the equalized signals and the reference constellation. However, de- 
pending on the constellation used during transmission, the PLL can 
have an A/-fold symmetric ambiguity, and thus it has limited CFO 
acquisition range; e.g., \ fk\ < 1/8 for 4Q AM signals. Moreover, 
the PLL require a long convergence time. To solve these problems, 
several methods have been proposed |3|, |5|, |6|, llll that allow 
for blind estimation of the CFO and symbols using the second-order 
cyclo-stationary statistics of the over-sampled received signal. Blind 
CFO estimation has also been studied in the context of orthogonal 
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frequency-division multiplexing (OFDM) systems, where the CFO 
destroys the orthogonality between the carriers (see f4| and the ref- 
erences therein). 

In a spatially distributed multiuser antenna system where data 
are transmitted simultaneously from multiple antennas, the received 
signal contains multiple CFOs, one for each transmit antenna. A 
PLL does not work in this case as there is no single frequency to 
lock onto. The literature on estimation of multiple CFOs is rather 
sparse. In |8|, multiple CFOs were estimated by using pilot sym- 
bols that were uncorrected among the different users. To account 
for multiple offsets, [10 J proposed that multiple nodes transmit the 
same copy of the data with an artificial delay at each node. The 
resulting system was modeled as a convolutive single-input/single- 
output (SISO) system with time-varying system response caused by 
the multiple CFOs. A minimum mean-square error (MMSE) deci- 
sion feedback equalizer was used to track and equalize the channel 
and to recover the input data. Training symbols were required in 
order to obtain a channel estimate, which was used to initialize the 
equalizer. 

In this paper, a new approach to blind CFO estimation and sym- 
bol recovery is proposed. The received base-band signal is over- 
sampled, and its polyphase components are used to formulate a vir- 
tual MIMO problem. By applying blind MIMO system estimation 
techniques, the system response can be estimated, and decoupled 
versions of the user symbols can be recovered, each one of which 
contains a distinct CFO. By applying a PLL, the CFO can be miti- 
gated and the transmitted symbols can be recovered. The estimated 
MIMO system response provides information about the CFOs that 
can be used to initialize the PLL, speed up its convergence, and avoid 
ambiguities usually linked with PLLs. 

2. SYSTEM MODEL 

We consider a distributed antenna system, where K users trans- 
mit simultaneously to a base station. Narrow-band transmission is 
assumed here, where the channel between any user and the base 
station is frequency non-selective. In addition, quasi-static fading 
is assumed, i.e., the channel gains remain fixed during the packet 
length. The continuous-time base-band received signal y(t) can be 
expressed as 

K 

y{t) = J2 '^*^*(* - r^y^^^-* + w{t) , (1) 

k = l 

where at represents the effect of channel fading between the k—th 
user and the base station and also phase offset; Tk is the delay asso- 
ciated with the path between the fc— th user and the base station; Fk 
is the frequency offset of the fc— th user and ^(t) represents noise; 



Xk{t) denotes the transmitted signal of user k: 

Xk{t) = ^Sk{i)p{t-iTs) , (2) 

i 

where Sfe(i) is the i— th symbol of user fc; Ts is the symbol period; 
and p(t) is a pulse function with support [0, Ts\. 

Our objective is to obtain an estimate of s(i) = 
[si(z), SK{iy\^ in the form 

s(i) = AP^s(i) , (3) 

where P is a column permutation matrix and A a constant diago- 
nal matrix. These are considered to be trivial ambiguities, and are 
typical in any blind problem. 

3. FORMATION OF THE MIMO PROBLEM 

The received signal y{t) is sampled at rate l/T = P/Ts, where 
the over-sampling factor P > /C" is an integer. In order to guar- 
antee that all the users' pulses overlap at the sampling times, the 
over-sampling period should satisfy: Ts/P > Tk,k = 1,...K. 
Or, in other words, the over-sampling factor P is upper bounded 
by Ts/ min{ri, ...,tk}- 

Let t = iTs + mT, m — 1, . . . , , P, denote the sampling times. 
The over-sampled received signal can be expressed as 

ym{i) = y{iTs+mT) 

K 

k=l 
K 

fe=i 

K 

= 5]a™,fc(s4«)e''"^"')+™(i+5), "i = l,...,P, (4) 

k-l 

where fu = FkTa is the normalized frequency offset between the 
fc— th user and the base station, and the typical element of the virtual 
MIMO channel matrix A is given by 

am,k = afce^'^*'> (^T, - r^) . (5) 

Define the following: y(i) = [yi{i), ...,yp(i)f \ 
A = {am,k}, a tall matrix of dimension P x K; 

s{i) = [si(i)e^'2"^i\...,sx(i)e^'2"-f^']^; and w(i) = 

[w{i + -p),...,w{i + -p)]"^. Then, ^ can be written in ma- 
trix form as 

y(j)=As(i)+w(i). (6) 

We could use the training based method of |8| to solve the 
MIMO system l|4j. That method assumes that the pilot symbols of 
different users are uncorrelated. The CFOs are obtained by search- 
ing for the location of a peak in the cross-correlation between the 
Discrete-Time Fourier Transform (DTFT) of a pilot sequence and 
that of the received signal. 

In the following we show how to estimate CFOs and recover the 
transmitted signals in a bind fashion, i.e., without the need for pilot 
symbols. The advantage of a blind approach is bandwidth efficiency 
since no bandwidth is wasted transmitting pilot symbols. 



4. BLIND CHANNEL ESTIMATION AND 
COMPENSATION OF THE CFOS 

Let us make the following assumptions. 

• Al) For each m. — 1, . . . , P, w,n{.) is a zero-mean Gaus- 
sian stationary random processes with variance a^, and is 
independent of the channel inputs. 

• A2) For each k, the sequence Sfc(i) is a zero 
mean with independent and identically distributed 
(i.i.d.) elements having nonzero kurtosis; i.e., 

jt^ = Cum[sk{i),sl{i),Sk{i),sl{i)] / 0. The sequences 
Sfe's are also mutually independent. 

• A3) The over-sampling factor P is no less than K. 

Under assumption (A2), it is easy to verify that the rotated input sig- 
nals Sfc(.) are also zero mean and i.i.d with nonzero kurtosis. Also, 
the Sk{i)'s are mutually independent for different fc's. Assumption 
(A3) guarantees that the virtual MIMO channel matrix A in l|6]( has 
full rank with probability one. If the delays of users are randomly 
distributed in the interval [0,Ts/P), then each row of the channel 
matrix can be viewed as having been drawn randomly from a contin- 
uous distribution so that the channel matrix has full rank with prob- 
ability one. 

One can apply any blind source separation algorithm (e.g., 
(T],|2lor(7l)toobtain 

A = APA . (7) 

Subsequently, using a least-squares equalizer we can obtain an 
estimate of the de-coupled signals s(i), within permutation and 
scalar ambiguities as 

i(i) = (A^A)"' A"y(i) = e^'^''''^-'^* | Ar'P^s(i) . (8) 

Without loss of generality we can assume that the transmitted signal 
has unit power. Then, on denoting by 6k the fc— th diagonal element 
of Arg{A}, the j— th separated input signal can be expressed as 

Ifc(i) =Sfc(j)e^(-*'=+'"^'''> . (9) 

In order to recover the transmitted signals, we still need to mit- 
igate the effect of CFO in each decoupled signal. This can be done 
via a PLL. By using the decoupled signals as inputs, and the constel- 
lation used in transmission as a reference, the PLL can effectively 
mitigate the CFO by minimizing the feedback error, which is calcu- 
lated based on the distance of the recovered signal and the closest 
valid constellation point. However, depending on the constellation 
used in transmission, there is a four-fold symmetric ambiguity for 
MQAM signals, or M -fold symmetric ambiguity for MPSK signals. 
For example, for 4QAM signals, and an initial CFO value of fk — 0, 
the effective tracking range for fk is < 1/8. Moreover, de- 
pending on the value of the CFO, the PLL generally needs a long 
convergence time, during which the input signals are not correctly 
recovered. 

Next we will show that by exploring the structure of the virtual 
channel matrix, we can obtain an estimate of the CFOs, which can 
then be used to initialize the PLL. By doing this, we can prevent 
the symmetric ambiguity problem and enlarge the effective tracking 
range of the PLL from \fk\ < 1/8 to |/fc| < 1/2. Also, the conver- 
gence time of the PLL can be greatly reduced. 



By taking the phase of the estimated channel matrix A, we ob- 
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where 0^ — Argjafe} + 9k, which accounts for both the phase of 
Ok and the estimated phase ambiguity in 

By applying linear fitting on the th column of 'i/ we obtain 
the least squares estimate of fj as 



(Ep.iP)(Ep=i*P 
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(11) 



We can write fj=fj+ where ej represents estimation error. 

On noting that the de-coupled signals Sj (i) in (|9} are shuffled in 
the same manner as the estimated CFOs in ilH . we can use the esti- 
mated CFOs to compensate for the effects of CFO in the decoupled 
signals ^ and thereby obtain estimates of the input signals as 



(12) 



Due to errors in the channel estimates, we can only compensate 
for most, but not all, of the CFO effects in (|9} and so we can write 



Sk{i) = Sk{i)e 



j(-«fe-2T6fci) 



(13) 



By subsequently applying a PLL to Sj{i), we can further mitigate 
the effects of the residual CFO e^. For 4QAM signals, as long as 
|efc| < 1/8, the residual CFO can be effectively removed by the PLL. 
The initial CFO estimator l lllb can prevent the symmetric ambiguity 
of the PLL, and can also greatly reduce the convergence time of the 
PLL. From dlOt , we can see that the CFO estimator will achieve full 
acquisition range for the normalized CFO. 

5. SIMULATION RESULTS 

In this section, we verify the validity of the proposed method via sim- 
ulations, under the following assumptions. The channel coefficients 
ak, k = 1, . . . ,K are zero-mean Gaussian random variables. The 
waveform p(-) is a Hamming window. The delays, Tk, k — 1, . . . , K 
are uniformly distributed in the range [0, Ts/P). The input signals 
are 4QAM signals. 

The blind source separation algorithm used here 
is the JADE method, which was downloaded from 
http://www.tsi.enst.fr/cardoso/guidesepsou.html 

First we show results for a two-user systems with /i — 
-0.1552, /a = 0.4335, ai = 0.3173 - 0.6483j, aa = 0.1625 + 
0.5867i, with SNR = 2QdB, and N = 1024. In Fig. Qwe show 
the polyphase outputs yi , 1/2 . Due to the mixing and the CFOs no 
obvious constellation is visible. In Fig. |2l we show the de-coupled 
signals s^, fc = 1, 2 right after JADE. Although still rotated by the 
CFOs, two signals s^, fc = 1, 2 are clearly separated. In Fig. [S] we 
show the recovered input signals Sk, k — 1,2, where we can see 
that after compensating for the effect of CFOs the constellations are 
recovered. 

Next we show estimation results averaged over 300 indepen- 
dent channel runs, and 20 Monte-Carlo runs for each channel. For 
each channel case, the coefficients ak,k = 1,2 were generated 
randomly, and the continuous CFOs where chosen randomly in the 



range [— 25^)- delays, Tk, k = 1, K where chosen uni- 
formly in the range of [0, Ts/P). The transmitted signal was 4Q AM. 

The performance of both the pilot-based method and the pro- 
posed method at different data lengths and with SNR set to 30dB 
is shown in Fig. [4] For the pilot-based method, each user trans- 
mits a pilot sequence of length 32, and the pilots are random se- 
quences uncorrelated between different users. Fig. |4] shows the 
mean-square error (MSE) for the CFO estimate ([tB for different 
values of the over-sampling factor P. The MSE is calculated usings 

^ Ef=i[(/fe - = w Ef=i[(A - Fk)Tsf. We can see that 
by increasing P we can improve the estimation accuracy. Fig. [5] 
shows the Bit Error Rate (BER) for different values of P. For both 
blind and training methods, the BER is calculated based on the re- 
covered signals after the PLL. As expected, the BER performance 
also improves by increasing P. The proposed method appears to 
work well even for short data length. 

Next we show the performance of both methods at various noise 
levels. We use packet length A'^ — 1,024. In Fig. |6]we show the 
MSE of the blind CFO estimator (TTJ as well as that of the training 
based method. We can see that by increasing P we improve estima- 
tion accuracy. In Fig. |7] we show the BER performance after the 
PLL for both blind and training based methods. We see that the pro- 
posed blind method has almost the same performance as the training 
based method for SNR values lower than 20dB, while the training 
based method can achieve better BER performance for higher values 
of SNR. 



6. CONCLUSION 

In this paper we have proposed a novel blind approach for identifica- 
tion of a distributed multiuser antenna system with multiple CFOs. 
By over-sampling of the received base-band signal, the MISO prob- 
lem is converted into a MIMO one. Blind MIMO system estimation 
then yields the system response, and MIMO input recovery yields 
the decoupled transmitted signals, each one containing a CFO. By 
exploring the structure of the MIMO system response we obtain a 
coarse estimate of the CFOs, which is then combined with a deci- 
sion feedback PLL to compensate for the CFOs in the decoupled 
transmitted signals. The proposed blind method has full acquisition 
range for normalized CFOs. 
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Fig. 1. Received mixing signal y 



Fig. 4. MSE of CFOs vs iV for K=2, with SNR=30dB, 4QAM 
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Fig. 5. BER vs N for K=2, with SNR=30dB, 4QA1VI 
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Fig. 7. BER vs SNR for K=2, 4QAM, T=1024 



